Nucleation of a non-critical phase in a fluid near a critical point 
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Abstract 

Phase diagrams of some globular proteins have a fluid-fluid transition as well as a fluid-crystal transition. Homo- 
geneous nucleation of the crystal from the fluid phase near the critical point of the fluid-fluid transition is examined. 
As the fluid-fluid critical point is approached, the number of molecules in the critical nucleus, the nucleus at the top 
of the free energy barrier to nucleation, is found to diverge as the isothermal compressibility. This divergence is due 
to a layer of the fluid phase of width equal to the fluid's correlation length which surrounds the core of the nucleus; 
the number of molecules in a crystalline environment in the nucleus does not diverge. The free energy barrier to 
nucleation remains finite but its derivative with respect to the chemical potential is equal to minus the number of 
molecules in the critical nucleus and so diverges. 

1 Introduction 

The phase behaviour of a number of globular proteins has been studied and in addition to the fluid-to-crystal transition 
there is a metastable fluid-fluid transition,^' ^ i.e., a transition between two fluid phases differing in density which ends 
at a critical point. This is analogous to the vapour-liquid transition of simple substances such as water. The transition 
is metastable: it exists within the fluid-crystal coexistence region. A fluid near the fluid-fluid critical point is anomalous, 
essentially because the correlation length is very large. Here we consider the free energy barrier Afl* to crystallisation 
in a fluid which is close to the critical point. Crystallisation starts by the nucleation of a microscopic crystallite, this 
then grows to form a crystal. The rate at which such microscopic crystallites form is proportional to exp(— Ar2*/fcT) 
because their formation is an activated process."^'* Here we study only the size and free energy of the critical nucleus, 
we leave consideration both of the dynamics of its formation and of its growth to later work. We also limit ourselves 
to temperatures above or equal to the critical temperature. Our main flnding is that the number of molecules in the 
critical nucleus varies as xt, the isothermal compressibility, near the critical point. As xt diverges at the critical point, 
so does the number of molecules in the nucleus. Within our theory the number of molecules in the critical nucleus is 
equal to minus the derivative of the free energy barrier with respect to the chemical potential and so near the critical 
!" ■ point the barrier to nucleation is a rapidly decreasing function of chemical potential. 
. ^ We focus on the universal aspects of critical nuclei near critical points.^ It is well-known that near a critical point 
■ fluids exhibit universal behaviour, behaviour which is determined solely by the universality class of the system,^ here 
^ ' that of the three-dimensional Ising model. Although in the centre of a critical nucleus the density will be far from the 
■ - - ' critical density, at the fringes of the nucleus far from the centre the density will be close to the bulk density of the fluid, 
and there in the fringes we find universal behaviour. We find that the fringe of the nucleus dominates the number of 
molecules in the nucleus but not the free energy. It does however dominate the derivative of the free energy with respect 
to the chemical potential as this is nothing but minus the number of molecules in the nucleus. In order to focus on 
the universal aspects we use a simple phenomenological theory and for simplicity we use a mean-fleld theory, although 
such theories have well-known deflciencies near critical points.^ We use a simple theory flrst used for droplets by Cahn 
and Hilliard.^ It is variously called Cahn-Hilliard theory, van der Waals-Cahn-Hilliard theory or the square gradient 
approximation.^' ^"^^ 

Numerical work on droplets near a critical point has been done by Talanquer and Oxtoby^-'^ using a theory similar 
to a Cahn-Hilliard theory. This followed pioneering computer simulations of a critical nucleus near a critical point by 
ten Wolde and Frenkel.^^ Talanquer and Oxtoby obtained nuclei with very large numbers of molecules near a critical 
point but did not perform the analytical analysis required to extract out the scaling near a critical point. Although 
their theory has an additional order parameter the scaling of the number of molecules near the critical point in their 
model is almost certainly the same as found here. See Refs. 13-15 for other recent theoretical work on nucleation 



near a metastable transition. Also, more than twenty years ago, Widom^' used Cahn-Hilliard theory to look at the 

(planar) interface between coexisting phases, where one of these phases was at a critical point. He envisaged not a 
single component system with fluid-fluid and fluid-crystal phase transitions but a binary mixture with liquid-liquid 
demixing and vapour-liquid phase transitions. However, the universal aspects are the same for both systems. Thus, our 
calculations for the radial density profile of the nucleus are analogous to those for density profile of the planar interface 
found by Widom. In addition, our results are relevant to the nucleation of the vapour phase near the critical point of a 
metastable liquid-liquid demixing transition. 

In the next section we derive analytic expressions, within a mean-field theory, for the density in the fringe of a nucleus 
in a near-critical fluid. Then we go on to derive an expression for the number of molecules, and hence the derivative of 
the free energy barrier near the critical point. The third and final section is a conclusion. 

2 Theory 

The critical nucleus is at the top of the barrier, it is at a maximum in the excess grand potential. The excess grand 

potential is the grand potential with a nucleus minus that without a nucleus. To start we require an expression for the 
excess grand potential of a nucleus Afl as a functional of the density function of the nucleus and then to find the density 
function which extremises this functional. For a crystalline nucleus, a rather complex functional is required, see for 
example the functional used to calculate the fluid-crystal interfacial tension for hard spheres. However, the universal 
effects we are looking for here occur at the fringe of the nucleus where the nucleus will have a density close to that of 
the fluid. There a much simpler functional suffices, such as those used to calculate fluid-fiuid interfacial tensions and 
the excess grand potentials of fluid droplets in fluid phases. 

For a spherically symmetric fluid droplet in a fluid, the density profile (/>(r) is a function only of the distance r from 
the centre of the droplet. The fluid phase has a density p and a chemical potential fi; <p is the density at a point whereas 
p is the bulk density of the fluid. Of course (f){r cx)) = p. The standard Cahn-Hilliard expression for the grand 
potential cost Aft of such a droplet, as a functional of its density profile (/'(r), is^-^-io 

dr, (1) 

where 

ALj{<t>) = m-f{p)-Kct>-p), (2) 

is the work required per unit volume to change the density from p to (j>, at a chemical potential p.. f{(p) is the bulk 
Hclmholtz free energy per unit volume of the fluid at a density (f). The second term in Eq. (1) within the brackets is 
the gradient term: the grand potential cost due to variations in space of the density. The gradient squared term is the 
lowest order term in a gradient expansion and so is adequate when the density is slowly varying. The coefficient, k of 
this term is taken to be a constant. See Refs. 3,9, 10 for its relation to the intermolccular potential. 

The functional Eq. (1) will be totally inadequate within the crystalline core of a nucleus but we will not require 
either the density function or the contribution to the grand potential of this core. Thus, we will use Eq. (1) as our 
functional and end up with expressions which are the sum of two terms, one term from the core for which our expression 
will be wrong but which we will not evaluate, and another term from the fringe, for which Eq. (1) is a reasonable 
(mean-field) approximation which we will evaluate. 

The critical nucleus is at the top of the free energy barrier and so is at a maximum of AQ. Thus, for the critical 
nucleus we may set the functional derivative of AQ with respect to the density profile (p{r) to zero, 

(f^)_..VV = 0. (3) 

The curvature of the density profile at a point is proportional to the derivative of the excess grand potential with respect 
to the density at that point. Using Eq. (2) for Au) in Eq. (3) we have 

/x(0) -p- 2kV^^ = 0, (4) 

where the p without an argument is the chemical potential in the fluid, and p{4>) is the chemical potential of a bulk 
fluid at a density (j). Once we have solved Eq. (4) we can insert the solution into Eq. (1) to obtain the excess grand 
potential of the critical nucleus, denoted by Afl* . 

The excess number of molecules in the critical nucleus. An*, is the number of molecules with the nucleus present 
minus the number without it. It may be obtained by integrating over the density profile, 



Afl 



- 



An* = I A(/)(r)dr, 



(5) 



where A(^(r) = (f){r) — p is the excess density at a point, the density at a point minus the bulk density. The derivative 
of ACl* with respect to the chemical potential of the fluid /x is, using Eq. (1), 



= - / A0(r)dr = -An*. (6) 

T J 

Thus, the derivative of Afl* with respect to the chemical potential is simply —An*. This result, that the derivative of 
the excess grand potential of the critical nucleus is minus the excess number of molecules in the nucleus, is often called 
the nucleation theorem. 



2.1 The fringe of the nucleus 

Equation (4) can be solved numerically if the chemical potential is known as a function of density. Here, wc would like 
to concentrate on the fringe of the critical nucleus. The fringe is the outermost part of the nucleus, where the density 
is near the density of the fluid, p. We define the fringe as being that part of the nucleus which is more than a distance 
Tc from the centre. The distance rc is such that A(j){r > rc)/p ^ 1. Fig. 3 is a schematic of the radial density profile of 
the nucleus. Now, if the fractional density difference A(j)/p <C 1 we can use a Taylor expansion for — p., 

where xt is the isothermal compressibility of the fluid (at a density p). The isothermal compressibility is defined as^^ 

Substituting Eq. (7) into Eq. (4) we have 

-^A(j)(r) - 2KV^A(j)(r) =0 A0 «: p. (9) 
P XT 

This is the Helmholtz equation, for Ac/), and the solution is a function of the Ornstein-Zernike form, 

A0(r) = p%xp(-r/0, (10) 
with ^ the correlation length of the fluid, given by 

e = 2kp\t. (11) 

To obtain Eq. (10) the boundary conditions A(j){r — > oo) and A0(rc) = p{a/rc) exp{—rc/^) were employed, a is 
a molecular length scale, a few nms for proteins. A^/p will, as required for Eq. (9), be small for r > Vc provided that 
Tc is a few times a or more. We do not need to specify Vc beyond saying that it must be at least a few times a. This 
implies a core a few molecules across, which is reasonable. From Eq. (10) we see that the width of the fringe is, as we 
might have expected, of the order of the correlation length ^ in the surrounding fluid. This width will thus diverge as 
the fluid approaches a critical point. 



2.2 Near a critical point 

We now consider a nucleus in a fluid which is near a critical point, either at the critical density but just above the critical 
temperature, or at the critical temperature but at a density near the critical density. In either case as the critical point 

is approached the isothermal compressibility, and hence (sec Eq. (11)) diverges. Within mean-field theory, along the 
critical isochore, i.e., with the density fixed at its value at the critical point pcp, Xt scales with temperature as^'^ 

XT~(T-Tej,)-i p = Pcp T>T^p, (12) 

where Tcp is the temperature of the critical point. The compressibility diverges as one over the temperature ciifFercnce 
to the critical point. An alternative path to the critical point is along the critical isotherm. We fix the temperature 
T = Tcp and vary the density. Along the critical isotherm xt varies as 



XT ~ (p - Pep)-' T = T, 



(13) 



The compressibility diverges as one over the square of the density difference to the critical point. 

Putting our solution for A(/)(r > rc), Eq. (10), into Eq. (5) wc find that when ^ is very large and so the fringe 
correspondingly large in volume, that the fringe dominates An*. Equation (5) is easily evaluated by substituting Eq. 
(10) for A(j){r > rc), 

/•r=rc 

An* = / A^(r)dr + Awpa^,'^ C > ''c 

Jr=0 
pr=rc 

= I A^(r)dr + Sttctk/O^Xt ^ » Tc, (14) 

ir=0 

where wc used the relation between ^ and xt, Eq. (11). The excess number of molecules An* is a sum of two terms. The 
first term, the integral in Eq. (14), comes from the core of the nucleus. It is of order {pc — p)r^ where pc is the density 
of the phase which is nucleating. It remains finite as the critical point is approached. However, the second term in Eq. 
(14), the contribution of the fringe to the number of molecules in the critical nucleus An*, diverges. As the critical 
point is approached the number of molecules in the nucleus tends to infinity, as suggested by Talanquer and Oxtoby.^^ 
It does so as the compressibility xt- This is the central result of this work, and is a general result independent of the 
nature of the phase which is nucleating. The fringe of the nucleus is at densities close to the fluid's density and far from 
that in the core of the nucleus. So, the divergence in An* is from a divergent number of molecules at densities near the 
bulk fluid density: although they are part of a nucleus of a crystalline phase they themselves are in a fluid environment. 
The number of molecules which are in a crystalline environment docs not diverge as the critical point is approached. 
Thus, within mean-field theory, along the critical isochore the number of molecules diverges as (T — Tcp)~^ and along 
the critical isotherm as {p — pcp)~'^ ■ 

Now we consider the contribution of the fringe to the grand potential of the critical nucleus. We require the Taylor 
expansion of Aw about its value at the bulk density, 

A.(« = i^ + .... (15) 

2 P^XT 



The first nonzero term in the Taylor expansion of Aw is the quadratic term as the first density derivative of Aw is zero 
at p. Now, splitting the integration in Eq. (1) at Tc, and using the quadratic approximation, Eq. (15), for Aw for r > Tc, 
we have 

■i(^ + K(VA^)^ 



Af2: 



/ ' [Aw + K (V(A)^l dr + /" 

Jr=0 '- Jr 



2 P^XT 



" ■ dr, (16) 



where, for r > rc, we also substituted A<p for <p in the gradient term. After substituting Eq. (10) for Ac/) in the second 
integral, we see that the contributions to the integral from both the terms in the integrand decreases as a function 
of radial distance r. Taking into account the factor of in dr the integrand varies as exp(— 2r/^) times a quadratic 
polynomial in 1/r: it is a monotonically decreasing function of r. The farther we go out into the fringe of the nucleus 
the smaller is the contribution to the excess grand potential, the barrier to nucleation. This means that the barrier to 
nucleation will be dominated by contributions from the core of the nucleus where the behaviour is non-universal and 
our approximations are invalid. Recall that our functional, Eq. (1), is a very poor approximation in the core and hence 
the first integral in Eq. (16) is not an accurate expression for the contribution of the core to the excess grand potential. 
We are thus unable to calculate the absolute excess grand potential of a critical nucleus near a critical point but we are 
of course able to calculate its derivative with respect to the chemical potential as this is just minus the excess number 
of molecules in the critical nucleus, which is dominated by the fringe. 



3 Conclusion 

Some globular proteins^' ^ have two transitions: a fluid-crystal transition, call it transition a, and a fluid-fluid transition, 

call it transition /3. Transition a is always strongly first order but transition f3 becomes continuoiis at the critical point. 
Here, we have studied homogeneous nucleation associated with transition a near a continuous transition /3. Systems near 
a continuous transition show universal (within a universality class) behaviour and are very susceptible to a perturbation 
— response functions such as the compressibility take large values. The nucleus associated with homogeneous nucleation 
of a transition a therefore perturbs the highly susceptible fluid near transition /?. Because the fluid is near this continuous 
transition, the response functions such as the isothermal compressibility xt, are large. The perturbation due to the 
nucleus is therefore large, it extends out to the large distance ^ and involves a number of molecules which scales with 
Xt- The free energy barrier to homogeneous nucleation of transition a is then a rapidly decreasing function of chemical 
potential. As suggested by ten Wolde and Prenkel,^^ the continuous transition (3 helps the dynamics of transition a. 
This is true in general, not just for crystallisation near a critical point. The nucleation of the vapour phase when a 



mixture of liquids boils near a metastable critical point of liquid-liquid demixing is another example. Note that there 

An* is negative and will tend towards — oo not +00. 

The theory of section 2 is a simple mean-field theory. However, our conclusion that An* ~ xt near the critical point 
may be valid beyond mean-field theory. If so then the divergence of An* along a particular path in the diagram, e.g., 
the critical isochore, will occiir with whatever is the correct critical exponent of xt along that path. For example, for 
nucleation near a critical point in the universality class of the Ising model in three dimensions, xt scales as (T — r^p)"^'^^ 
along the critical isochore not (T — Tcp)~^ as mean-field theory predicts, see Refs. 3,6. 

The starting point for most theoretical studies of nucleation is classical nucleation theory.'^''* Within classical 
nucleation theory the critical nucleus is modeled as a small crystallite, assumed to have bulk properties, separated from 
the surrounding metastable fluid by an interface, assumed to have the same surface tension as that between the coexisting 
bulk phases and to be thin, no more than a couple of molecules thick. For a nucleus in a fluid far from a fluid-fluid 
critical point, this assumption of a thin interface with the bulk surface tension is not unreasonable. The correlation 
lengths in the fluid and crystalline phases will both be small and so will the width of the interface. Then if the width of 
the interface is small with respect to the diameter of the nucleus then approximating its surface tension by that of the 
planar interface between coexisting phases is reasonable. However, for a nucleus in a near-critical fluid, the correlation 
length in the surrounding fluid ^ may be very large. So, for a crystal nucleus which is only a few lattice spacings across 
the interface between it and the surrounding fluid may be much thicker than the diameter of the crystalline part of the 
nucleus. Then this interface will be very different from that assumed within classical nucleation theory, and so classical 
nucleation theory will yield only a poor estimate of the barrier to nucleation. 
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